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Abstract 

A search for single top quark production is performed in e^p collisions at HERA. The 
search exploits data corresponding to an integrated luminosity of 118.3 pb -1 . A model for 
the anomalous production of top quarks in a flavour changing neutral current process in- 
volving a twy coupling is investigated. Decays of top quarks into a b quark and a W boson 
are considered in the leptonic and the hadronic decay channels of the W. Both a cut-based 
analysis and a multivariate likelihood analysis are performed to discriminate anomalous 
top quark production from Standard Model background processes. In the leptonic channel, 
5 events are found while 1.31 ± 0.22 events are expected from the Standard Model back- 
ground. In the hadronic channel, no excess above the expectation for Standard Model pro- 
cesses is found. These observations lead to a cross section a(ep — * etX) = 0.29 i q 14 
at yfs = 319 GeV. Alternatively, assuming that the observed events are due to a statistical 
fluctuation, upper limits of 0.55 pb on the anomalous top production cross section and of 
0.27 on the twy coupling K tui are established at the 95% confidence level. 
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1 Introduction 



In ep collisions at the HERA collider, the production of single top quarks is kinematically pos- 
sible due to the large centre-of-mass energy, which is well above the top production threshold. 
In the Standard Model, the dominant process for single top production at HERA is the charged 
current reaction e + p — > vtbX (e~p — > vtbX). This process has a tiny cross section of less than 
1 fb EEl and thus Standard Model top production is negligible. However, in several extensions 
of the Standard Model, the top quark is predicted to undergo flavour changing neutral current 
(FCNC) interactions, which could lead to a sizeable top production cross section. FCNC inter- 
actions are present in models which contain an extended Higgs sector [3|, Supersymmetry 01, 
dynamical breaking of the electroweak symmetry [5j or an additional symmetry 0. An obser- 
vation of top quarks at HERA would thus be a clear indication of physics beyond the Standard 
Model. 

The HI Collaboration has reported 000 the observation of events with energetic isolated 
electrons and muons together with missing transverse momentum in the positron-proton data 
collected between 1994 and 2000. The dominant Standard Model source is the production of 
real W bosons. However, some of these events have a hadronic final state with large transverse 
momentum, which is atypical of W production. These outstanding events may indicate a pro- 
duction mechanism involving processes beyond the Standard Model. One such mechanism is 
the production of top quarks which predominantly decay into a b quark and a W boson. The 
lepton and the missing transverse momentum would then be associated with a leptonic decay 
of the W boson (W — > tv), while the observed high P T hadronic final state would be produced 
by the fragmentation of the b quark. 

In this paper we present a search for anomalous single top production using the HI detector. 
The analysis uses the data collected between 1994 and 2000 corresponding to an integrated 
luminosity of 118.3 pb -1 . The search covers the leptonic decay channel (W — > iv) and the 
hadronic decay channel (W — > qq') of the W boson that emerges from the top quark decay. 
Both a cut-based analysis and a multivariate likelihood analysis are performed to select top 
quark candidates. Another search for single top production at HERA was recently published by 
the ZEUS Collaboration 0. 



2 Phenomenology of FCNC Top Production 

A FCNC vertex involving the direct coupling of the top quark to a light quark (u or c) and a 
gauge boson would lead to single top production, as illustrated in figure [T] The most general 
effective Lagrangian, proposed in ifTOl . which describes FCNC top quark interactions involving 
electroweak bosons is: 



r FCNC 
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U =u,c 

9 



2cos# 
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where = (i/2) [7^, 7^], 9 W is the Weinberg angle, e and g are the couplings to the gauge 
groups with U{\) and SU{2) symmetries, respectively, t\j is the electric charge of up-type 
quarks, and Z^ are the fields of the photon and the Z boson and A is the scale up to which 
the effective theory is assumed to be valid. The tensors A^ v and Z^ v are the field strength tensors 
of the photon and Z boson fields. By convention A is set equal to the top mass in the following. 
Because gauge invariance is required, only magnetic operators allow FCNC couplings of the top 
quark to a photon and an up-type quark U — u,c, denoted by n tUl , while the non-vanishing Z 
mass allows both the magnetic coupling K tuz and the vector and axial vector couplings to a Z 
boson and an up-type quark, denoted by v tuz and a tuz . 

In ep-collisions, due to the large Z mass, the contribution of the Z boson and the 7 — Z 
interference are highly suppressed. Single top production is thus dominated by the t-channel 
exchange of a photon. Therefore Z-exchange is neglected in this analysis and only the n t u 1 
couplings are considered. 

In the Standard Model FCNC processes can only arise via higher order corrections and are 
strongly suppressed. The possibility of anomalous single top production at HERA was first 
investigated in H, where a model was built in which K tqi oc m 2 q and therefore only the tcy 
coupling was relevant. In order to cover a large range of possible underlying theories, the 
model described by CQ) allows both couplings n tcl and K tul to be present and independent of 
each other. The sensitivity of HERA is much larger to the coupling K tul than to n tcJ since the 
w-quark density in the proton is much larger than the c-quark density at the high Bjorken x 
values needed to produce top quarks. In the analysis presented in this paper only the coupling 
K tul to the w-quark is considered. 

Compared with the production of top quarks in a 7-u fusion process, the corresponding 
charge conjugate anti-top quark production in a j-u fusion process involving sea quarks is 
suppressed by a factor of ~80. It is not considered in this analysis. 

The simulation of the anomalous single top signal relies on an event generator (ANOTOP), 
which uses the leading order (LO) matrix elements of the complete e+q — > e+t — > e+b+W — ► 
e + b + f + f' process as obtained from the CompHEP [111 program. The BASES/SPRING O 
package is used to perform the numerical integration of the amplitudes and to generate events 
according to the resulting differential cross section. The parton shower approach lfP31 . which 
relies on the leading logarithmic DGLAP [fT4l evolution equations, is used to simulate QCD 
corrections in the initial and final states. The MRST LO parton densities are used for the 
proton lfT5l . The parton densities are evaluated at the top mass scale. The nominal top mass 
is set to 175 GeV. A variation of the top mass by +5 GeV (—5 GeV) induces a cross section 
variation of —20% (+25%). The cross section calculation for anomalous single top production 
has recently been improved by including next-to-leading order (NLO) QCD corrections |[T6l . 
The uncertainty related to the choice of the renormalisation and factorisation scales is reduced 
to about 5%. These NLO QCD corrections increase the cross section by about 17% and are 
taken into account as an overall correction factor to the LO calculation for the results derived in 
this analysis. 
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3 Standard Model Background Processes 



Signatures of single top production are searched for in the leptonic and hadronic decay channels 
of the W boson that emerges from the top quark decay. The relevant final state topology for 
the leptonic channel is an isolated lepton with high transverse momentum, at least one jet and 
missing transverse momentum. For the hadronic channel, the signature is three or more jets 
with high transverse momenta. The Standard Model processes that produce events with similar 
topologies and thus constitute the background for the present analysis are outlined below. Due 
to its small cross section, Standard Model top production is not considered in this analysis. 

The Standard Model processes that produce background to the leptonic channel were inves- 
tigated in detail in the context of the isolated electron and muon analyses published by the HI 
Collaboration in [ 8 1 and are only briefly described here. The main contribution is the produc- 
tion of W bosons with subsequent leptonic decay of the W. The production of the electroweak 
vector bosons W ± is modelled using the EPVEC ifPTI generator. The NLO QCD corrections 
to W production lfT8ll are taken into account by weighting the events as a function of the ra- 
pidity and transverse momentum of the W boson |[T9"ll . Other processes can contribute to the 
investigated final state through misidentification of photons or hadrons as leptons or through 
fake missing transverse momentum due to measurement fluctuations. Processes with a genuine 
lepton but possible fake missing transverse momentum are lepton pair production and neutral 
current (NC) deep inelastic scattering. The contribution from lepton pair production, dominated 
by two-photon processes where one of the two produced leptons is not detected, is calculated 
with the GRAPE [ 20 1 generator. The background contribution from NC deep inelastic scat- 
tering is estimated using the RAPGAP lETll generator. In charged current (CC) deep inelastic 
scattering, the missing transverse momentum is genuine, but a hadron or a photon from the final 
state may be falsely identified as a lepton. This background contribution is calculated using the 
DJANGO C2I program. 

For the hadronic channel the production of multi-jet events in photoproduction and NC deep 
inelastic scattering are the most important Standard Model backgrounds. The RAPGAP ET1 
generator is used to model multi-jet production in NC deep inelastic scattering for virtualities 
of the exchanged photon Q 2 > 4 GeV 2 . Multi-jet events with photon virtualities Q 2 < 4 GeV 2 
are generated with the PYTHIA Monte Carlo program ll2"31 . Both generators rely on first order 
QCD matrix elements and use leading-log parton showers and string fragmentation lfT3l . Both 
light and heavy quark flavours are generated. The GRV LO (GRV-G LO) parton densities [1241 
in the proton (photon) are used. The production of W bosons and their hadronic decay are 
simulated using the EPVEC generator. This contribution is negligible in the present analysis. 

All generated events are passed through the full GEANT ll25ll based simulation of the HI 
apparatus and are reconstructed using the same program chain as for the data. 



4 Experimental Conditions 

The analysis is based on e ± p collisions recorded by the HI experiment between 1994 and 2000. 
At HERA electrons or positrons with an energy E e of 27.6 GeV collide with protons at an 
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energy of 920 GeV, giving a centre-of-mass energy of ^fs = 319 GeV. Up to 1997 the proton 
energy was 820 GeV, giving y/s = 301 GeV. The data correspond to an integrated luminosity 
of 37.0 pb 1 in e+p scattering at- s /s = 301 GeV <X 30 i = 37.0 pbr 1 ), together with 13.6 pb 1 
in e~p scattering and 67.7 pb -1 in e + p scattering at ^fs = 319 GeV (£319 = 81.3 pb -1 ). 

A detailed description of the HI detector can be found in [26|. Only those components 
essential for this analysis are briefly described here. The right handed Cartesian coordinate 
system used in the following has its origin at the nominal primary ep interaction vertex. The 
proton direction defines the z axis. The polar angle, 9, and transverse momenta, P?, are defined 
with respect to this axis. The region 9 < 90° is referred to as the "forward" region. The 
pseudorapidity is defined as 77 = — ln(tan|). 

The inner tracking system contains the central (25° < 9 < 155°) and forward (7° < 9 < 
25°) drift chambers. It is used to measure the trajectories of charged particles and to determine 
the position of the interaction vertex. Particle transverse momenta are determined from the 
curvature of the trajectories in a solenoidal magnetic field of 1.15 Tesla. 

Hadronic final state particles as well as electrons and photons are absorbed in the highly 
segmented liquid argon calorimeter [1271 (4° < 9 < 154°), which is 5 to 8 hadronic in- 
teraction lengths deep depending on the polar angle. It includes an electromagnetic section 
which is 20 to 30 radiation lengths deep. Electromagnetic shower energies are measured with 
a precision of a(E)/E = 12%/,/E/GeV © 1%, hadronic shower energies with a(E)/E = 
50%/ y 1 Ej GeV © 2%, as determined in test beam measurements [28|. In the backward re- 
gion (153° < 9 < 178°), the liquid argon calorimeter was complemented by a lead-scintillator 
backward electromagnetic calorimeter (BEMC) before 1995 and by a lead-scintillating fibre 
spaghetti calorimeter (SPACAL) J29| since 1995. 

The calorimeter is contained within a superconducting coil and an iron return yoke, instru- 
mented with streamer tubes, which is used as a muon detector and covers the range 4° < 9 < 
171°. Tracks of penetrating particles, such as muons, are reconstructed from their hit pattern in 
the streamer tubes and are detected with an efficiency of above 90%. The instrumented iron also 
serves as a backing calorimeter to measure the energies of hadrons that are not fully absorbed 
in the liquid argon calorimeter. 

In the forward direction, muons are also detected in the forward muon system, a set of drift 
chambers covering the range 3° < 9 < 17°. This detector measures the muon track momenta 
from their curvature in a magnetic field provided by a toroidal iron magnet. 

The trigger conditions for interactions leading to high transverse energy in the final state, as 
expected for top quark production, are mainly based on liquid argon calorimeter signals. Events 
in the leptonic channel are triggered by their calorimetric missing transverse momentum. The 
trigger efficiency is 50% (85%) for events with a missing transverse momentum above 12 GeV 
(25 GeV). Events containing an electron with an energy of at least 10 GeV are triggered via 
the energy deposition in the electromagnetic calorimeter with an efficiency larger than 95%. 
Events with muons may also be triggered by a set of triggers based on signals consistent with a 
minimum ionising particle in the muon system in coincidence with tracks in the inner tracking 
system. In the hadronic channel, the triggering of events with three or more high P T jets is based 
on the scalar sum of the transverse energy deposited in the liquid argon calorimeter. For events 
containing three jets with transverse momenta above 25 GeV, 20 GeV and 15 GeV respectively, 
the trigger efficiency is close to 100%. 
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5 Search for Single Top Production in the Leptonic Channel 



The decay cascade t — > bW — > 6£z/ yields events with a lepton, missing momentum and a jet. 
The search for top quarks starts with the selection of events containing a high Pt lepton and 
missing transverse momentum. In these events kinematic reconstruction of potential top quark 
decays is performed yielding a preselected data sample which serves as the basis for the final 
selection of top quark candidates. To discriminate single top production from Standard Model 
background processes, observables characteristic of top quark decays are used in a cut-based 
top selection and later in a multivariate likelihood analysis. 

5.1 Events with isolated leptons and missing transverse momentum 

The search for top quark decays in the leptonic channel is based on the selection of events with 
high Pt leptons and missing transverse momentum described in [8|. Further details on this 
analysis can be found in OU1I3T1I . The selection yields a sample of candidates for leptonic W 
boson decays. The following variables are used to characterise the events. 

• P T : the transverse momentum of the lepton. Electron transverse momenta are calculated 
using calorimetric information together with vertex information from the trackers. Muon 
transverse momenta are measured from the curvature of the charged track detected in the 
central tracker or in the forward muon detector. 

• 9 e : the polar angle of the lepton. 

• Charge of the lepton: The charge is measured from the track associated with the lepton. 
It is considered to be determined if the signed curvature of the track is different from zero 
with a measurement accuracy of better than two standard deviations. For less accurate 
measurements the lepton charge is labelled as "undefined". 

• P™ JSS : the total missing transverse momentum reconstructed from all observed final state 
particles. 

• P*: the transverse momentum of the hadronic final state. The hadronic final state, de- 
noted by X, is reconstructed by combining energy deposits in the calorimeter and charged 
tracks as described in 02"ll . It does not include the energy deposited by any identified lep- 
tons in the event. 

The main selection criteria are the requirement that there be an electron or muon with high 
transverse momentum P T > 10 GeV in the polar angle range 5° < 9 e < 140° and large missing 
transverse momentum p™ ss > 12 GeV as a signature of the undetected neutrino from the W 
decay. The lepton is required to be isolated from neighbouring tracks or jets. The distance 
Dtrack m pseudorapidity-azimuth space of the closest track from the lepton is required 
to be > 0.5 and that of the nearest jet, Dj et , to be > 1.0. In the muon channel, an additional 
cut on the transverse momentum of the hadronic final state, P* > 12 GeV, is applied. Further 
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selection criteria are applied to suppress processes where leptons are faked or missing transverse 
momentum is induced by measurement fluctuations, as discussed in detail in [ 8 1. 

In the full e ± p data sample, 19 events |0 are selected, compared with a Standard Model 
prediction of 14.5 ± 2.0, where the latter is dominated by W production (10.7 ± 1.8). One of 
the 19 events was observed in e~p collisions. 

5.2 Kinematic reconstruction of the top quark decay 

In order to calculate the kinematics of a top quark decay, a reliable reconstruction of both the b 
quark and the neutrino from the W decay is necessary. 

Reconstruction of the b quark 

In most cases the b quark manifests itself as a single high P T jet in the detector. However, gluon 
radiation from the b quark can lead to final states with more than one jet. Therefore the b quark 
momentum is reconstructed as the sum of the momenta of all jets found in the event. Jets are 
identified using an inclusive k T algorithm ll3*3l with a minimum jet transverse momentum of 
4 GeV. The sum of all jets gives a better approximation to the b quark momentum than the full 
hadronic final state X, since it is less sensitive to particles originating from the proton remnant. 

Reconstruction of the neutrino 

The transverse momentum vector of the neutrino corresponds to the vector of the total missing 
transverse momentum: P£ = P™ tss . Concerning the longitudinal momentum of the neutrino, 
two cases are treated separately. 

• Tagged events: the scattered beam electron is detected in one of the calorimeters. This is 
expected to be the case for 30% of single top events. In the electron channel the scattered 
beam electron is assumed to be the one with the lower transverse momentum, which 
according to the simulation is the correct choice in 95% of top events. From energy and 
longitudinal momentum balance we obtain 

(E - P z ) u = 2E e -{E- p z ) le P tons -(E- P z f . (2) 

• Untagged events: the scattered beam electron is lost in the beam pipe and therefore its 
longitudinal momentum is unknown. In this case a constraint is applied on the invariant 
mass of the lepton and the neutrino from the W decay: 

M lv = ^Pf + PZ + 2P t P v « y/2P e P v = M W = 80.42 GeV [34j, (3) 

where Pi and P v denote the four-vectors of the lepton and the neutrino, respectively. 
The constraint on the W mass generally yields two possible solutions for (E — P z ) v . If 
two solutions for the neutrino kinematics exist, one solution corresponds to a backward 
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neutrino with 6 e < 9 V , where 9 V is the neutrino polar angle, while the other solution cor- 
responds to a forward neutrino with 9 e > 9 U . The solution which is most likely according 
to the model for anomalous top production is chosen. It is found that in top quark de- 
cays where the charged lepton is observed at small polar angles (9g < 18°) the backward 
neutrino solution is favoured, while for large lepton polar angles (fig, > 40°) the forward 
neutrino solution is favoured 051 . In the intermediate region of the polar angle neither 
solution is favoured. Therefore, for 18° < 9g < 40°, the solution is selected that yields 
an invariant mass Mg vh of the system consisting of the lepton, neutrino and b candidate 
closest to the nominal top mass of m t = 175 GeV. 

This neutrino reconstruction method is discussed in detail in ll35l . It has a reconstruction 
efficiency of 99% (95%) for simulated top events in the electron (muon) channel. The widths 
from Gaussian fits of the top reconstructed mass M&,& distributions obtained for simulated top 
decays are 13 (18) GeV for the electron (muon) channel. 

The kinematics of all 19 isolated lepton events allow the reconstruction of a neutrino ac- 
cording to the above procedure. In figure El the invariant mass M ivh of these events is plotted 
against the lepton-neutrino transverse mass, defined as: 

Mr =yJ{Pb + PtY ~ (Pt + Pt) 2 , 

where Pj. and Pj, are the transverse momentum vectors of the lepton and the neutrino respec- 
tively. For each event the measured lepton charge is also indicated. Several events are situated 
at large masses close to the top mass and have transverse masses Mjf compatible with 
there being a W boson in the event. 

5.3 Selection of top quark decays in the leptonic channel 

Top preselection 

In addition to the kinematic reconstruction of the top quark decay, the lepton charge is also 
exploited. The decay chain t — ► bW + — > b£ + vi produces only positively charged leptons. The 
production of anti-top quarks, which would yield negatively charged leptons, is strongly sup- 
pressed, as mentioned in section|2l To reduce the contribution from processes other than FCNC 
top quark production, negatively charged leptons are therefore rejected. The "top preselection" 
thus consists of the following three steps: 

• selection of isolated lepton events with missing transverse momentum; 

• neutrino reconstruction; 

• rejection of leptons with negative charge. 

The preselected top sample contains 9 electron events and 6 muon events compared with an 
expectation from Standard Model processes of 8.40 ± 1.06 and 1.88 ± 0.32, respectively. For 
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negative lepton charges, 4 events are found while 3.48±0.53 are expected from Standard Model 
sources. The event yields, Standard Model predictions and selection efficiencies after each step 
of the top preselection are summarised in table [T] 

Candidates for single top event production are searched for in the preselected top sample. 
Both the cut-based top selection described below and the multivariate likelihood analysis de- 
scribed in section|7]exploit kinematic observables that are characteristic for top quark decays. 

Observables for top quark decays in the leptonic channel 

The following three observables are chosen to discriminate single top production from Standard 
Model W production and other Standard Model processes. 

1. P|.: the transverse momentum of the b candidate. 

2. M luh : the invariant mass of the system consisting of the lepton, neutrino and b candidate, 
which corresponds to the top quark candidate mass. 

3. 9yy : the W decay angle - defined as the angle between the charged lepton momentum in 
the rest frame of the W boson and the W direction in the rest frame of the top quark. 

Distributions of the observables P|> , and cos 9^ are shown in figure |3] for the top 
preselection. Also shown are the signal distributions of simulated top events with an arbitrary 
normalisation. An excess of events at the highest values of P^ is visible in both the electron and 
the muon channels. The masses for some of the electron and muon events are compatible 
with the top quark mass. 

Cut-based top selection in the leptonic channel 

Starting from the top preselected sample, the cuts that are used to select top candidates are 
Pj. > 30 GeV and M ivh > 140 GeV. In this cut-based selection, no restriction is imposed on 
cos Oyy , since it does not yield an efficient separation of single top quark from Standard Model 
W production, while the contribution from other Standard Model processes is already reduced 
to a negligible level by the and M^f, cuts. 

In this cut-based analysis, 3 electron events and 2 muon events are selected as top quark 
candidates in the full e ± p data sample. Some properties of these events are presented in table|2| 
The background expectation from Standard Model processes is 0.65 ± 0.10 events for the elec- 
tron channel and 0.66 ± 0.12 for the muon channel, as summarised in table[TJ The efficiency for 
simulated top events is 36% (38%) for the electron (muon) channel, taking into account the top 
decays where the W boson decays via W — > r — > e(/x). Combining both channels, there are 5 
top quark candidates in the data for an expectation of 1.31 ± 0.22 from background processes. 

The systematic uncertainties on the Standard Model prediction that are relevant for the lep- 
tonic channel are described in [8|. They are dominated by the uncertainty of 15% on the NLO 
cross section calculation for Standard Model W production 1 18 ]. 
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6 Search for Single Top Production in the Hadronic Channel 



A search for single top production is also performed in the hadronic channel. The decay cas- 
cade t — > bW — > bqq' yields events with at least three jets with high transverse momenta. The 
main Standard Model background is the QCD production of high Pt jets in photoproduction 
and neutral current deep inelastic scattering. First, a high statistics sample of multi-jet events 
is compared with simulations of Standard Model multi-jet production to check that the back- 
ground is well understood. The search for top quarks is then performed in a multi-jet sample 
restricted to large transverse momenta ("top preselection"). Further details of the analysis can 
be found in ll35l . 



6.1 Multi-jet events 

Since only LO Monte Carlo simulations with leading-log parton showers are used to model 
QCD multi-jet production, these simulations may give only approximate descriptions of the 
shape and normalisation of kinematic distributions. Therefore, as a first step, the agreement in 
shape of the data and the simulations is studied using multi-jet events. The LO background 
simulations are normalized to the observed number of events in the data, to also account for the 
higher order QCD effects. 

Jets are reconstructed using the inclusive kx algorithm based on calorimetric energy de- 
posits combined with well measured tracks. To ensure a reliable measurement, only jets in the 
pseudorapidity range —0.5 < rf et < 2.5 and with transverse momenta P T et > 4 GeV are con- 
sidered. In order to remove electrons which are misidentifled as jets, each jet is required to have 
either an electromagnetic energy fraction of less than 90% or a jet size larger than 0. 1 , where the 
jet size is defined to be the energy weighted average distance in the r]-(p plane of the particles 
composing the jet from the jet axis. 

Events with at least three jets with P T etl > 25 GeV, P T et2 > 20 GeV and P T et3 > 15 GeV 
are used to study the agreement between the data and the simulations. The multi-jet data sample 
contains 1472 events. In figure distributions of various kinematic quantities are shown for 
these events and compared with the Standard Model predictions and a simulation of single top 
production. The transverse momenta of the three highest Pt jets, the total hadronic transverse 
energy (E^), the di-jet invariant mass closest to the W mass (M]J e c ° mfc ) and the invariant mass 
of all jets (Mj ets ) show good agreement with the Standard Model simulations in shape. 

The appropriate overall normalisation factors for the PYTHIA and RAPGAP simulations 
are determined using two complementary subsamples of the multi-jet events. One subsam- 
ple contains events where no electron is identified (low Q 2 sample), where the prediction is 
dominated by the PYTHIA simulation. The other subsample contains events with an identified 
electron (high Q 2 sample), in which case the prediction is dominated by the RAPGAP simu- 
lation. A normalisation factor of 1.29 is applied to the event yield predicted by PYTHIA for 
Q 2 < 4 GeV 2 and a factor of 1.40 is applied to RAPGAP for Q 2 > 4 GeV 2 . 
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6.2 Selection of top quark decays in the hadronic channel 

Top preselection 

The search for top quarks in the hadronic channel is performed in a sample which is further 
restricted to the high transverse momentum region defined by P^ etl > 40 GeV, P^ et2 > 30 GeV 
and P^ et3 > 15 GeV. Since top quarks typically deposit a large amount of transverse energy in 
the detector, a cut on the total hadronic transverse energy of E^ 1 > 110 GeV is also applied. In 
addition, one of the jet pairings must yield an invariant mass between 65 GeV and 95 GeV, cor- 
responding to a window around the nominal W mass with a width of twice the mass resolution 
obtained for hadronic W decays. These selection criteria are referred to as the top preselec- 
tion in the following. In the data 92 events are selected. After application of the normalisation 
factors, the expectation for Standard Model processes is 92.4 ± 16.6. The good agreement be- 
tween the data and the prediction in the top preselection indicates that the normalisation factors 
obtained for the multi-jet sample are also valid at high transverse momenta. 

Observables for top quark decays in the hadronic channel 

The observables used for the discrimination of the top signal from the QCD background are 
chosen in analogy to the leptonic channel (section 1531) . The jet among the three highest P T jets 
that is not used to form the W mass is assigned to the b quark ("6 candidate"). A study using 
simulated top events shows that this hypothesis correctly identifies the b quark jet in 70% of the 
events. The three characteristic observables used are: 

1. Pj,\ the transverse momentum of the b candidate; 

2. Mj ets : the mass of the top quark - reconstructed as the invariant mass of all jets in the 
event. The width of a Gaussian fit to the mass distribution obtained for simulated top 
decays is 14 GeV; 

3. '■ the W decay angle - defined as the angle in the W rest frame between the lower 
P T jet of the two jets associated to the W decay and the W direction in the top quark rest 
frame. The helicity structure of the decay implies that the lower P T jet corresponds to the 
q from the W decay in most of the top events. 

Distributions of these three observables are shown in figure |5] for the top preselection and are 
compared with the Standard Model processes and with the simulated top signal. Good agree- 
ment between the data and the Standard Model simulations is seen for all three distributions. 
No sign of an excess compatible with single top production is visible. 

Cut-based top selection in the hadronic channel 

In the hadronic channel harsher cuts need to be applied than in the leptonic channel to enhance 
the top signal. To select top quark candidates following the top preselection, the transverse 
momentum of the b candidate has to fulfill P? > 40 GeV and the invariant mass of all jets has 
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to be reconstructed in a window around the top mass given by 150 < Mj ets < 210 GeV. In 
addition, cos 8^ > —0.75 is required. For this selection the efficiency estimated with simulated 
top events is 30%. The number of candidate events selected is 18, compared with 20.2 ± 3.6 
events expected from Standard Model processes. In figure |5] the reconstructed mass Mj ets of 
the selected top quark candidates is shown. The observed data events can be accounted for by 
Standard Model processes. 

The main experimental systematic uncertainties on the expected number of events for the top 
selection are due to the uncertainty in the absolute hadronic energy calibration of the calorimeter 
(±4%) and the measurement of the polar angles of the jets (±20 mrad), leading to a total 
experimental uncertainty of 11%. Systematic effects due to the uncertainty in the luminosity 
measurement and the trigger inefficiencies are negligible. The normalisation of the Standard 
Model simulation is taken from the data and its uncertainty of 10% thus corresponds to the 
statistical uncertainty on the observed number of data events. An additional uncertainty of 10% 
accounts for differences between the shapes of the kinematic distributions in the data and the 
simulations. All uncertainties presented above are added in quadrature. The total systematic 
uncertainty on the expected number of events in the hadronic channel amounts to 18%. 



7 Multivariate Likelihood Analysis 

In addition to the cut-based top quark selections presented in sections 15.31 and 16.21 a multi- 
variate analysis is performed as an alternative approach to the search for single top production. 
The observables used to discriminate the top signal from the Standard Model background are 
combined to form a single discriminator based on the relative likelihood approach as defined 
in Poll. 

In this framework a set of i observables V = {Vj} with the corresponding densities p^ 9nal 
and p^ ack 9 round ^ calculated from Monte Carlo samples for the signal and background respec- 
tively, are used for each event to calculate a discriminator: 

'Tjsignal 

i 

Here C(V) denotes Gaussian correction factors used to correct for correlations between the 
variables Vi as explained in detail in Poll . The discriminator V(V) is an approximation to the 
likelihood that an event is part of the signal rather than the background. 

The variables used in this analysis are V = {P|., M eub , cos 6^} for the electron and muon 
channels and V = {P?, Mj ets , cos 6^} for the hadronic channel. The distributions of the dis- 
criminator variables for the signal and background processes according to the simulations and 
the distributions of the data events are shown in figure |6] for the electron, muon and hadronic 
channels. By definition, the top signal populates the region close to V = 1, while the Standard 
Model background populates the region close to V = 0. In both the electron and the muon 
channels, there are two populations in the data: one class of events that are more Standard 
Model-like and another class of events that are more top-like. The five top candidates selected 
in the cut-based analysis in section 15.31 correspond to the events with the highest likelihood 
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V > 0.7. In the hadronic channel the distributions for the data events and for the QCD expecta- 
tion are in good agreement. All candidates selected by the cut-based analysis in section l6"2l are 
situated at V > 0.35, but there is no evidence for any enhancement over the QCD background 
in this region. A comparison between the likelihood discriminator analysis and the results of the 
cut-based analyses (sections l5.3l and l6"2l) for each channel is presented in table[3]by applying a 
cut on the discriminator which yields the same signal efficiency as the cut-based selection. 

The probabilities for the Standard Model to fluctuate to discriminator distributions at least as 
unlikely as those observed in the data are evaluated using Monte Carlo experiments, following 
the approach proposed in 071 . They are 0.3% for the combined electron and muon channels, 
40.1% for the hadronic channel and 2.6% for the combination of all channels. 

The discriminator distributions are used to quantify the possible top signal contribution in 
the data using a maximum-likelihood fit. A likelihood function L is introduced as the product of 
Poisson probabilities of observing n k data events in each bin k of the discriminator distribution: 

n bin n/; 

L = T[e-^^-, (4) 

where fi k = B k + S k is the sum of the signal contribution S k and expected background B k 
in bin k. The total top signal normalisation, S = J2k £>k> is fitted as a free parameter, while 
the normalisation of the background is fixed to the Standard Model prediction. The value of 
S which best matches the data can be obtained by maximising the likelihood function L, or 
correspondingly by minimising the negative log-likelihood function —2 InL. By using the 
factor 2, the log-likelihood function corresponds to a x 2 function in the Gaussian limit. 

The log-likelihood functions after subtraction of the minimum values are shown in figure 13 
as functions of the single top cross section for the combined electron and muon channels, the 
hadronic channel and the combination of all channels. The conversion of the fitted signal nor- 
malisation S to a single top cross section at yfs = 319 GeV is done for each channel by folding 
in the corresponding efficiency e top , the top and W branching ratio product B t ^bw • Bw-+ff an d 
taking into account the integrated luminosities £ 30 i and £319: 

a(y/s = 319 GeV) = — 1 -— . (5) 

£to P ■ Bt-^bw " Bw-¥ff 0.70 • £301 + £319 

Here, the factor 0.70 is the ratio of the cross sections at ^Js = 301 GeV and 319 GeV fl6l . The 
branching ratio for t — > bW is assumed to be B t ~+bw = 100%, in accordance with |38|. The 
combination of the different channels is performed by adding the log-likelihood functions of 
the single channels. In order to propagate the systematic uncertainties related to the measure- 
ment through the signal fitting procedure, each observable affected by a systematic uncertainty 
is smeared according to a Gaussian distribution with a width corresponding to the size of the 
uncertainty. The full analysis is then repeated for a large number of values of each smeared 
observable. The r.m.s of the resulting distribution of shifts in the fitted cross section is taken as 
the corresponding systematic uncertainty. The uncertainties on the cross sections are added in 
quadrature for each error source. In order to include these uncertainties in the log-likelihood 
function, the function is approximated by a half-parabola on each side of the minimum (Gaus- 
sian approximation) and the width of the parabola is increased according to the total systematic 
uncertainty. As can be seen in figure|7J the impact of the systematic uncertainties is most impor- 
tant for the hadronic channel, while it is negligible compared with the statistical uncertainties 
in the leptonic channels. 
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8 Cross Sections and Limits 



The log-likelihood function for the combination of the electron and muon channels in figure 13 
yields a single top cross section of 0.411qa9 pb at ^fs = 319 GeV. This cross section is different 
from zero by more than two standard deviations and reflects the observation of the five top quark 
candidate events presented in section |5J For the hadronic channel the log-likelihood function in 
figure[7]yields a single top cross section of 0. 04^^23 pb at = 319 GeV. This cross section is 
consistent with no top signal, in accordance with the results of the cut-based analysis reported 
in section |6"21 Taking into account the statistical and systematic uncertainties, the results from 
the hadronic channel and the combined electron and muon channels are compatible at the 1.1 a 
level. The combination of all three channels yields a single top cross section of 0.29+QJ4 pb at 
yfs = 319 GeV. The addition of a contribution from anomalous single top production to the 
Standard Model background provides a good description of the data. As discussed in sectional 
the cross section for anomalous top production at HERA is approximately proportional to the 
anomalous magnetic coupling squared, a(ep — > e + 1 + X) oc k^ u . The value obtained for K tul 
is 0.20l°;°jj. 

In view of the small number of top candidates an upper limit on the single top production 
cross section is calculated. It can be directly obtained from the log-likelihood functions pre- 
sented in figure 13 A one-sided exclusion limit at the 95% confidence level corresponds to an 
increase of 2.69 units in — 2AlnL. The resulting upper bound on the single top cross section at 
y/s = 319 GeV for the combination of all analysed channels is 

a(ep -> e + t + X, = 319 GeV) < 0.55 pb (95% CL). 

The bound on the top cross section is translated into an upper limit on the anomalous twy 
coupling of: 

\K tul \ < 0.27 (95% CL). 

The limits obtained separately for the combined electron and muon channels and for the hadronic 
channel only are given in table |4j As a cross check, other statistical methods have also been 
used to derive the exclusion limits, for example a Bayesian approach with a flat prior ll34l or 
likelihood-based approaches 07113*91 . The results are consistent within 15%. 

The present limit on n tul is consistent with the result obtained by the ZEUS collabora- 
tion [9 1 in the framework of the NLO QCD calculation ifTBI for anomalous single top produc- 
tion: \ntwi\ < 0-17 (95% CL). The limits on the top quark anomalous couplings obtained by 
the HERA experiments can be compared with the limits obtained from the search for single 
top production at LEP iRUll and from the analysis of radiative top decays by the CDF collab- 
oration [38 1 . Figure E] represents the current status of the constraints on K tul and v tu z- The 
limit on the anomalous coupling K tul obtained in the present analysis significantly improves 
the CDF and LEP upper bounds if the vector coupling v tu z is not too large. The error band on 
the HI limit represents the uncertainty induced by a variation of the nominal top quark mass 
of m t = 175 GeV by ±5 GeV in the analysis. Other theoretical errors are neglected. The HI 
results on single top production are not in contradiction with the limits set by other experiments. 
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9 Summary 



A search for single top production is performed using the data sample collected by the HI 
experiment between 1994 and 2000, corresponding to a total luminosity of 118.3 pb -1 . This 
search is motivated by the previous observation of events containing an isolated lepton, miss- 
ing transverse momentum and large hadronic transverse momentum, a topology typical of the 
semileptonic decay of the top quark. 

In a cut-based analysis, 5 events are selected as top quark candidate decays in the leptonic 
channel. The prediction for Standard Model processes is 1.31 ± 0.22 events. The analysis 
of multi-jet production at high P T , corresponding to a search for single top production in the 
hadronic channel, shows good agreement with the expectation for Standard Model processes 
within the uncertainties. 

In order to extract the top quark production cross section, a multivariate likelihood anal- 
ysis is performed in addition to the cut-based analyses. The top signal contribution in each 
channel is determined in a maximum-likelihood fit to the likelihood discriminator distribu- 
tions. The results from the hadronic channel do not rule out a single top interpretation of the 
candidates observed in the electron and muon channels. For the combination of the electron, 
muon and hadronic channels a cross section for single top production of a = 0.29^o;|i pb at 
y/s = 319 GeV is obtained. This result is not in contradiction with limits obtained by other 
experiments. The addition of a contribution from a model of anomalous single top production 
yields a better description of the data than is obtained with the Standard Model alone. 

Assuming that the small number of top candidates are the result of a statistical fluctuation, 
exclusion limits for the single top cross section of a < 0.55 pb at y/s = 319 GeV and for the 
anomalous twy coupling of \n tui \ < 0.27 are also derived at the 95 % confidence level. The 
HERA bounds extend into a region of parameter space so far not covered by experiments at 
LEP and the TeVatron. 
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Electron Channel 


Data 


Standard Model 


W only 


Top efficiency 


Top 

preselection 


Isolated Lepton + P^ iss 


11 


11.53 ±1.49 


8.17 ±1.35 


47% 


v reconstruction 


11 


11.06 ± 1.43 


7.93 ± 1.31 


46% 


Cut on lepton charge 


9 


8.40 ± 1.06 


5.72 ±0.95 


45% 


Top cut-based selection 


3 


0.65 ±0.10 


0.57 ±0.10 


36% 


Muon Channel 


Data 


Standard Model 


W only 


Top efficiency 


Top 

preselection 


Isolated Lepton + P^ lss 


8 


2.96 ±0.50 


2.54 ±0.49 


46% 


v reconstruction 


8 


2.70 ±0.46 


2.38 ±0.46 


44% 


Cut on lepton charge 


6 


1.88 ±0.32 


1.67 ±0.32 


43% 


Top cut-based selection 


2 


0.66 ±0.12 


0.59 ±0.12 


38% 



Table 1 : Observed and predicted numbers of events for the three steps in the top preselection 
and for the cut-based top selection in the leptonic channel. The "W only" column gives the 
prediction from Standard Model W production alone. The numbers are presented for the full 
e ± p data sample corresponding to an integrated luminosity of 118.3 pb -1 . 
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Run 


Event 


Lepton 
type 


Charge 


P» 
(GeV) 


M eub 

M w -constraint solution 
(GeV) 


M% 
(GeV) 


248207 


32134 


e 


+ (15(7) 


43 


155+ 7 , 


63+2 


252020 


30485 


e 


+ (40a) 


47 


168i| 


5H1 


268338 


70014 


e 


+ (5.1(7) 


48 


160+J 


88+J 


186729 


702 


n 


+ (4.0a) 


72 




43lg 


266336 


4126 


n 


+ (26a) 


55 


172l? 


69^ 



Table 2: Kinematics and lepton charges of the five top quark candidates in the leptonic channel. 
In one event (Run 252020 Event 30485) the scattered electron is detected and the "tagged" mass 
solution can be obtained for the top quark mass 1681}} GeV. The mass of the lep ton-neutrino 
system of this event is measured to be Mg v = 79t\l GeV. 





Cut-based analysis 






Efficiency 


Data 


SM 


Data 


SM 


Electron Channel 


3 


0.65 ±0.10 


3 


0.67 ±0.13 


0.72 


36% 


Muon Channel 


2 


0.66 ±0.12 


2 


0.62 ±0.12 


0.40 


38% 


Hadronic Channel 


18 


20.2 ±3.6 


20 


17.5 ±3.2 


0.58 


30% 



Table 3: Observed and predicted numbers of events in the cut-based top selection, compared 
with the selection using the single cut on the likelihood discriminator (V > V min ) which yields 
the same efficiency for the top signal. The numbers are presented for the full e ± p data sample 
corresponding to an integrated luminosity of 118.3 pb -1 . 





(7 




Electron+Muon Channel 


< 0.90 pb 


< 0.35 


Hadronic Channel 


< 0.48 pb 


< 0.25 


All Channels 


< 0.55 pb 


< 0.27 



Table 4: Exclusion limits at the 95% confidence level for the single top cross-section at ^Js = 
319 GeV and for the anomalous tu^j coupling. 
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Figure 1 : Anomalous single top production via a flavour changing neutral current interaction at 
HERA, with subsequent decays t — > bW + and W + — > //'. 
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Figure 2: The invariant mass plotted against the lepton-neutrino transverse mass 
for the isolated electron (a) and muon (b) events after the neutrino reconstruction. The data 
events (points) are compared with simulated top events produced via FCNC interactions (small 
points, arbitrary normalisation) and events from Standard Model W production (open circles) 
corresponding to 50 times the integrated luminosity of the data. For each data event the lepton 
charge is indicated, if it is determined with a measurement accuracy of better than two standard 
deviations. 
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ELECTRON CHANNEL MUON CHANNEL 




Figure 3: Distributions of the observables P|>, M tl , b and cos 9^ for the top preselection in 
the electron channel (left) and the muon channel (right). In each figure the arrows indicate 
the measured values for the data events: 9 events in the electron channel and 6 events in the 
muon channel. The solid histogram corresponds to the total Standard Model expectation. The 
hatched histogram represents the contribution from Standard Model W production. The dashed 
histogram shows the distribution for simulated top events with an arbitrary normalisation. 
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HADRONIC CHANNEL - MULTI-JET EVENTS 
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Figure 4: Distributions of the multi-jet events with P 

r>jet3 
T 



,jetl 
T 



100 150 200 250 



M jets (^V) 



> 25 GeV, P 3 T et2 > 20 GeV and 



> 15 GeV The jet transverse momenta are shown for the three highest P T jets, as well 
as the total hadronic transverse energy, E^, the invariant di-jet mass closest to the W mass, 
M^f rnb , and the invariant mass of all jets, Mj ets . The data (symbols) are compared with the 
Standard Model simulations (histograms) after application of the normalisation factors of 1 .29 
(1.40) for PYTHIA (RAPGAP). The expectation from Standard Model processes is dominated 
by low Q 2 multi-jet production. The DIS contribution for Q 2 > 4 GeV 2 is shown as a hatched 
histogram. The error band represents a systematic uncertainty of 18% on the total Standard 
Model prediction. The expected shape of the top signal is shown as the dashed histogram with 
an arbitrary normalisation. 
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HADRONIC CHANNEL - TOP PRESELECTION 
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Figure 5: Distributions of the observables P^, cos 6^ and Mj ets for the top preselection in the 
hadronic channel. The data (symbols) are compared with the Standard Model simulations (his- 
tograms) after application of the normalisation factors of 1.29 (1.40) for PYTHIA (RAPGAP). 
The expectation from Standard Model processes is dominated by low Q 2 multi-jet production. 
The DIS contribution for Q 2 > 4 GeV 2 is shown as a hatched histogram. The error band repre- 
sents the systematic uncertainty of 18% on the total Standard Model prediction. The expected 
shape of the top signal is shown as the dashed histogram with an arbitrary normalisation. The 
Mj ets distribution is also shown after the full cut-based top selection (lower right). 
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Figure 6: Distributions of the discriminator V for the data candidates (symbols) and the Stan- 
dard Model expectation (histograms) for the top preselection in the electron, the muon and 
the hadronic channels. The solid histogram represents the Standard Model background. The 
dashed histogram shows the distribution for simulated top events with an arbitrary normalisa- 
tion. For the hadronic channel, the top signal is also shown normalised according to the cross 
section derived in the combined electron and muon channels (hatched histogram) and added to 
the Standard Model histogram. 
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Figure 7: The zero-suppressed log-likelihood function — 2AlnL as functions of the single top 
cross section at ^Js = 319 GeV for the combined electron and muon channels, the hadronic 
channel and the combination of all channels. For the latter, the one-sided (upper) exclusion 
limit at the 95% confidence level on the single top cross section is marked by the dashed-dotted 
line. 
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Figure 8: Exclusion limits at the 95% confidence level on the anomalous tqj magnetic coupling 
Kt ul and the vector coupling v tu z obtained at the TeVatron (CDF experiment Il38l0 . LEP (L3 
experiment is shown, which currently gives the best limit of the LEP experiments |40|) and 
HERA (HI and ZEUS experiments). The anomalous couplings to the charm quark are neglected 
K tey = v tc z = 0. The error band on the HI limit shows the uncertainty on the coupling n tul 
induced by a variation of the nominal top quark mass by ±5 GeV. 
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